This line only printed with preprint option 

Fission properties of the BCPM functional 

Samuel A. Giulianj*] and Luis M. RobledJll 

Departamento de Fisica Teorica, Universidad Autonoma de Madrid, E-28049 Madrid, Spain 

Abstract 

We explore the properties of the Barcelona Catania Paris Madrid (BCPM) energy density func- 
tional concerning fission dynamics. Potential energy surfaces as well as collective inertias relevant 
in the fission process are computed for several nuclei where experimental data exists. Inner and 
outer barrier heights as well as fission isomer excitation energies are reproduced quite well in all 
the cases. The spontaneous fission half lives isFare also computed using the standard semiclasical 
approach and the results are compared with the experimental data. A reasonable agreement with 
experiment is found over a range of 27 orders of magnitude but the theoretical predictions suffer 
from large uncertainties associated to the values of the parameters entering the spontaneous fission 
half life formula. The impact that increasing the pairing correlations strengths has in the sponta- 
neous fission half lives is analyzed and found to be large in all the nuclei considered. Given the 
satisfactory description of the trend of fission properties with mass number we explore the fission 
properties of the even-even uranium isotope chain from 226 U to 282 U. Very large half lives are found 
when getting close to neutron number N=170. 
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I. INTRODUCTION 



Fission is an unique physical phenomenon taking place in heavy atomic nuclei and in- 
volving the evolution of the system from its ground state to scission going through different 
intrinsic shapes that cover a wide range of different intrinsic deformation parameters jl-3|. 
Fission properties depend upon the competition between the surface energy term coming 
from the strong nuclear interaction and the Coulomb repulsion and therefore fission proper- 
ties are often used as constraints to refine the parameters of effective nuclear interactions. 
Typ,ca. examples are the Gogny Q D1S parametritic that was fine tuned to reproduce 
the fission barrier of 240 Pu [5]. More recently it has been used with Skyrme like interactions 
to define the UNEDEF1 parametrization 6|, [7|. The gross features of fission can be under- 
stood from a mean field perspective involving the Hartree Fock Bogoliubov (HFB) theory 
{9), 0] and therefore it not surprising the large amount of studies devoted to this subject 
with Skyrme like interactions |t 10 1 , Gogny ones 



field 



7fl or based on the relativistic mean 
20]. Fission observables also depend on the inertia of the system to the relevant 



collective degrees of freedom and therefore they are sensitive to pairing correlations. As 

a consequence, fission is a good testing ground to test both the theories and interactions 

commonly used in nuclear structure. In addition, the theoretical understanding of fission 

is relevant to other areas like energy production with nuclear reactors, radioactive waste 

degradation or the nucleosynthesis of heavy elements in the explosive galactic environments 

and following the r-process. Last but not least, a better understanding of fission could open 

the door to a better estimation of magic numbers and hence extra stability of super-heavy 

nuclei beyond Z=114. In this paper we explore the ability of a newly proposed energy den- 

r~ n 

sity functional (EDF) denoted as Barcelona-Catania-Paris-Madrid (BCPM) 123 to describe 



fission. The BCPM 



23] is a recent parametrization of the BCP EDF 21 



22, 24 



29| for 



nuclear structure calculations. Its free parameters have been adjusted to binding energies 
of even-even nuclei all over the nuclide chart, including deformed ones. Instead of the more 
traditional approaches, like Skyrme or Gogny like [jj], where some central potential form is 
guessed (contact, gaussian, Yukawa, etc) and used afterwards to fit nuclear matter proper- 
ties and/or the nuclear matter equations of state EoS (both symmetric and neutron), in the 
BCPM functional we start from a microscopic EoS that is fitted by means of a low order 
polynomial in the density. That polynomial fit is translated to a finite nuclei EDF just by 
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replacing the nuclear matter density by the density of the finite nucleus. This procedure is 
inspired by the local density approximation LDA and is common practice in condense matter 
physics in practical applications of the Kohn-Sham theory. The EDF is complemented with 
a finite range surface term, a contact spin-orbit interaction of the same form as in popular 
Skyrme or Gogny forces, the Coulomb interaction and finally, a density dependent zero range 
pairing interaction [27] with strengths fitted to reproduce Gogny's neutron matter pairing 
gap. The parameters of the functional (essentially those of the finite range surface term plus 
some freedom in the polynomial fit to fine tune the binding energy per nucleon) are fitted 
to reproduce binding energies of the 518 even-even nuclei of the 2003 mass table evaluation 
of Audi and Wapstra. The properties of the interaction concerning quadrupole, octupole 
and fission dynamics have also been explored j^]. As shown below, the BCPM functional 
gives reasonable results for fission observables like spontaneous fission half lives, fission iso- 
mer excitation energies and inner and outer barrier heights. We have also shown that those 
results could be improved by slightly modifying the amount of pairing correlations, either 
by modifying the pairing strengths or by going beyond the mean field approximation to 
restore particle number symmetry broken by the HFB procedure. As a consequence of the 
satisfactory performance of BCPM in describing fission, we have explored fission properties 
of the uranium isotopic chain from proton drip line to the neutron drip line. 



II. METHODS 

The BCP energy density functional |24j] is made of a bulk part which is inspired by fully 
microscopic and realistic calculations of symmetric and neutron matter equations of state 



25 



26] . The two equations of state (symmetric and neutron matter) given as a a function 



of the nuclear density are parametrized using a fitting procedure in terms of low order 
polynomials of the densities. To account for finite size effects related to the surface energy, 
a phenomenological finite range interaction of gaussian type is included. In addition, the 
Coulomb interaction and the spin-orbit term are taken exactly as in the Skyrme or Gogny 
forces. To deal with open-shell nuclei we include in the BCP functionals a zero range density- 
dependent pairing interaction fitted to reproduce the nuclear matter gaps obtained with the 
Gogny force (27|. The calculations in finite nuclei are carried out with a modification of the 
HFBaxial computer code |s8 1 . 
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To describe fission we have resorted to the traditional procedure based on the mean field 
approach with pairing correlations: the Hartree- Fock- Bogoliubov (HFB) with constraints. 
As constraining operators we have used mainly the axially symmetric quadrupole moment 
operator Q 2 o = z 2 — ^(x 2 + y 2 ) although some exploratory calculations have also been 
performed with the octupole moment operator Q30 arid the necking operator Qn{zq) — 
exp( — (z — Zq) 2 /Cq). Axial symmetry has been preserved in the calculations because of the 
high computational cost involved in the release of this restriction. We are aware of the 
relevance of triaxiality specially in the height of the inner fission barrier but its effect is 
merely quantitative and to a much lesser extent, qualitative. On the other hand, reflection 
symmetry is allowed to break at any stage of the calculation leading to octupole deformed 
shapes and asymmetric fission. As a consequence of the breaking of the parity symmetry 
we are forced to constraint the center of mass to the origin as to prevent spurious solutions 
corresponding to a translation of the nucleus as a whole. The quasiparticle creation and 
annihilation operators of the HFB theory are expanded in a harmonic oscillator basis (HO) 
preserving axial symmetry and containing HO states with J z quantum numbers up to 35/2 
and up to 25 quanta in the z direction. The two lengths of the HO basis have been optimized 
in a few nuclei and the b± and b z lengths obtained for each value of the quadrupole moment 



used 
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in nearby nuclei. The minimization process uses the gradient method as described in 



31] . The most evident advantage of this method is its handling of constraints, which 



allows a larger number of them to be treated at once. 

Dynamical parameters like quadrupole collective masses and moments of inertia are com- 
puted both in the so-called Adiabatic Time Dependent Mean Field Framework (ATDHFB) 
and the generator coordinate method (GCM) (see 35|, [36] for specific details ) . In both 
theories the evolution of the inertia parameters involve the inversion of huge matrices re- 
lated to the stability matrix of the HFB theory and therefore we have recurred to the usual 



approximations of neglecting the off diagonal 



stability matrix (cranking approximation 
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(corresponding to the interaction) terms in the 



34]). Zero point energy corrections to the HFB 



energy eo(Q2o) are also considered in the ATDHFB and GCM approaches. In addition, the 
rotational energy correction computed following the phenomenological prescription of Ref 



32] is also subtracted. This correction is very important to the shape of the potential energy 



as its value increases with deformation and can reach several MeV for large deformations. 
The spontaneous fission half life is computed with the standard WKB formalism of quan- 
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turn mechanics. In the WKB formalism the t s { is given (in seconds) by the formula 

t rf = 2.861(r 21 (l + exp(2S)) 
The action S along theQ2o constrained path is given by 

rb 

S= / dQ 20 pB(Q w )(V(Q w ) - (E GS + E 0) ). 

where the integration limits correspond to the classical turning points for the energy Eqs + 
E . For the collective quadrupole inertia B(Q 2 o) we have used both the ATDHFB and 
the GCM expressions. The results obtained with the two different theories can differ in 
several orders of magnitude as the ATDHFB mass is known to be a factor in between 1.5 
and 2 larger than the GCM mass. The potential energy V(Q2o) is given by the HFB mean 
field energy with two quantities subtracted: the zero point energy £0(^20) associated to the 
quadrupole collective motion and the rotational energy correction given by the rotational 



approximation with the Yoccoz moment of inertia (see [32| for further details). Finally, an 
additional parameter Eq is introduced in the action and is added to the ground state energy 
Eqs • It is meant to represent the quantal ground state energy obtained after considering 
quantal fluctuations in the quadrupole degree of freedom around the ground state. This 
quantity could be estimated to be half of the square root of the curvature around the 
minimum divided by the collective inertia but it is often taken as a free parameter or kept 
fixed to some reasonable value. We have followed the later approach with E = 1.0 MeV 
and estimated the impact of considering a larger value by repeating the calculations with 
Eq = 1.5 MeV. The spontaneous fission half life obtained in this way is subject to several 
uncertainties that can lead to differences of several orders of magnitude. The uncertainties 
are: 1) the height of the inner fission barrier gets reduced when triaxiallity is considered in 
the mean field calculation. The amount of reduction is typically of a couple of MeV ( but it 



can show some isotopic variation, see [llj for a recent account in the actinide region). 2) The 
value of Eq can also make a difference in tsF specially for long lived isotopes where the fission 
barrier is wide. The reason is that the value of Eq alters the classical turning points. 3) 
The values of the correlation energy corrections to the HFB energy included in V(Q2o)- The 
values are computed under certain assumptions and approximations and a better estimation 
of their values can lead to some changes to V(Q2q). 4) The approximations involved in the 
evaluation of the collective masses can lead to differences of the order of 40 or 50 % in their 
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values 5) Finally, the amount of pairing correlations is also an important ingredient both 
in the evaluation of the zero point energy corrections as in the evaluation of the collective 
inertia. As shown below, changes of a few percent in the pairing strength values can lead to 
changes in the theoretical estimation of the half lives in the range of five to seven orders of 
magnitude. 

On the other hand, the experimental values of the parameters defining the potential 
energy of the fission process, namely the inner and outer barrier heights and the excitation 
energy of the fission isomer seem to be more robust quantities to compare with as they are 
not as sensitive to pairing correlations as the other parameters. However, these experimental 
quantities are obtained under some model dependent assumptions that can mask the physical 
meaning of some of the parameters. Therefore, although we have compared our values with 
the experimental ones, we prefer to compare with spontaneous fission life times and more 
precisely with the evolution with Z and N of those quantities. We will also study the impact 
of changing various quantities entering the WKB formula. 



III. RESULTS 



A. Comparison with other interactions 

Before comparing with the experimental data a comparison with the HFB results obtained 
with the two Gogny interactions, namely, D1S and DIM, is in order. The Gogny D1S 
interaction has been used in a thorough study of heavy nuclei properties, including fission, 
in Ref jll| and it has proved to reproduce quite nicely most of the properties analyzed. 
On the other hand, the fission properties of DIM have not been analyzed in detail but its 
good properties in many other aspects of nuclear structure like binding energies, radii and 
collective mode properties make it a good candidate to compare with. As we have already 
made comparison with D1S concerning fission properties [23] of actinides ( 240 Pu) and super- 
heavies ( 262 Sg), we will just explore another actinide: the nucleus 234 U. In the lowest panel of 
Fig[T]we compare the HFB energy as a function of Q20 f° r the three functionals in the nucleus 
234 U. We observe how the shape of the three curves look rather similar up Q20 = 60 b, apart 
from a constant shift of a few MeV, and from there one the DIM interaction produces a less 
pronounced decline than D1S and the BCPM results are in between but closer to DIM than 
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Figure 1: Comparison of HFB mean field quantities as a function of the mass quadrupole mo- 
ment Q20 obtained with three different interactions. The BCPM EDF and the D1S and DIM 
parametrizations of the Gogny force. In the lowest panel, the HFB energy is given. In the middle 
panel the particle-particle correlation energy E pp = — Tt(Ak) is plotted for protons (dashed lines) 
and neutrons (full lines) for the three different kinds of calculations. In the panel above the oc- 
tupole and hexadecapole deformation parameters are given. Finally, in the upper-most panel the 
ATDHFB collective inertia is given. 



to D1S. In the next panel a quantity that can be used as a measure of pairing correlations, 
the particle-particle correlation energy given by E pp = — Tt(Ak) , is depicted as a function 
of the mass quadrupole moment for both protons (dashed lines) and neutrons (full lines). 
Again the shape of the curves for protons and neutrons look rather similar in the three cases 
but there is again a shift in E pp such that BCPM yields lower values than D1S and DIM. 
Those lower values correspond to less intense pairing correlations in the BCPM case being 
the quenching of proton pairing correlations more intense than the neutron ones. It is worth 
noticing that the pairing correlations in the proton channel for BCPM are significantly less 
intense that for the Gogny interactions. In the next panel, the octupole and hexadecapole 
moments as a function of Q20 are a l so given in the three cases. The results in the three cases 
are very similar, to the extent of appearing as a single curve for many Q20 values. Finally, 
in the upper-most panel the collective inertia in the ATDHFB approximation is plotted. 
As in the case of the HFB and particle-particle energies the shape of the curves are quite 
similar but the BCPM inertia is up to a factor of three larger than the inertias obtained 
with the Gogny forces. The large value of the BCPM inertia is a direct consequence of 
the quenched pairing correlations: less pairing correlations imply a lower gap and therefore 
smaller two quasiparticle energies. As the two quasiparticle energies enter the collective 
inertia in the denominator, quenched pairing correlations imply enlarged collective inertias. 
On the other hand, the DIM inertia is around 15% smaller than the D1S one, consistent 
again with the quenched pairing correlations in DIM as compared to the D1S ones. The 
GCM inertias, not depicted, look rather similar in shape to the ATDHFB ones but are a 
factor of 0.5-0.6 smaller. If the t$F half life is computed with the WKB formula for the three 
cases with the GCM inertias, the values t SF = 1.5 x 10 40 s, 6.6 x 10 34 s and 5.7 x 10 23 s are 
obtained for BCPM, DIM and D1S, respectively. The large differences observed of up to 7 
orders of magnitude can be attributed partly to the difference in the HFB curve but mostly 
to the different values of the collective inertias. The previous values have been obtained 
without taking into consideration the reduction of the inner barrier height as a consequence 
of triaxiallity. Also, increasing the value of the ^parameter reduces the half lives by 6, 2 
and 4 orders of magnitude respectively. Given the large variability of the half lives with the 
different parameters entering the WKB formula a direct comparison with the experimental 
data is meaningless and only comparisons along a series of nuclei, all of them computed with 
the same conditions, can lead to meaningful conclusions regarding fission properties. 
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B. Varying pairing strengths 



In the BCPM functional the pairing interaction is taken as a density dependent contact 
pairing interaction with strength parameters fixed to reproduce the neutron matter pairing 
gap of the Gogny force 27] . We have shown in the previous subsection that the particle- 



particle correlation energy, a quantity related to the amount of pairing correlations, was 
much smaller for BCPM than for the Gogny forces leading to much larger collective inertias. 
It is therefore reasonable to investigate the behavior of fission properties as a function of 
the pairing strength. To this end, a parameter r] has been introduced as a multiplication 
factor in front of the pairing gap field Ajt/. For the sake of simplicity we have considered 
an unique parameter for both protons and neutrons although different parameters will give 
more freedom to adjust to experimental data. The outcome of the calculations with rj values 
of 1.05 and 1.10 for the nucleus 234 U are presented in Fig[2j We observe in the lowest panel 
how increasing the pairing strength by 10 % (rj — 1.10) we gain of the order of 1 MeV in 
binding energy and the gain is even larger on top of the inner barrier. The net effect is 
to slightly decrease the inner barrier eight (Bj) whereas the other parameters, namely the 
outer barrier height Bjj and the fission isomer excitation energy remain more or less the 
same. The particle-particle correlation energies E pp for protons and neutrons increase with 
increasing r\ but the slope is larger for neutrons than for protons. The multipole moment 
values do not change at all in they are on top of each other for different 77 values. Finally, 
the impact on the collective inertia is clearly visible. Increasing the pairing strength by 5 % 
reduces the collective inertia by roughly 30 % whereas a 10 % increase leads to a reduction 
of a 50 %. The impact on the tsF is dramatic, decreasing by 11 orders of magnitude in going 
from rj = 1.0 to r\ — 1.05 and six additional orders of magnitude in going from r\ = 1.05 
to rj = 1.10. This result is a clear indication of the very important role played by pairing 
correlations in the description of fission. The result suggests that experimental fission data 
could be used to fine tune the pairing strength instead of more traditional approaches based 
on odd-even staggerings. From a theoretical perspective the result also points to the very 
important role that the correlations associated to particle number symmetry restoration 
should have in fission dynamics. Restoring particle number symmetry usually leads to 
larger pairing correlations than the ones present at the mean field level and therefore will 
have a tremendous impact on fission half lives. In this respect it is worth mentioning that 
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the dependence on density of BCPM is on integer powers of the density allowing the use 
of the regularization techniques suggested to solve some technical problems associated to 
;he evaluation of energy kernel overlaps typical of symmetry restoration theories (see Refs 



371 . |38|and references therein). 



C. Nuclei with known experimental data 

In order to validate the BCPM functional as a functional able to describe fission, we have 
performed calculations for those even-even nuclei where the spontaneous fission half life 
has been measured. We will also compare the parameters defining the theoretical potential 
energy surface, namely the inner and our barrier heights (Bi and Bjj) and the excitation 
energy of the fission isomer Ejj with available experimental data. It has to be mentioned 
that the experimental data for Bj and Bu is model dependent and therefore less reliable 
that the pure tsF data. Recently, a thorough comparison of theoretical results with the 



experimental data has been carried out in Ref. 
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. In table[T]the experimental and theoretical 



39 



40]. 



values for Bf, Bn and En are given for all nuclei where experimental data exists 
The theoretical values have been obtained by considering the HFB energy as a function of 
Q20 with the rotational energy correction (computed in the way described in the previous 
section) subtracted. The effect of the zero point energy correction e (Q2o) has not been 
included mainly because it is almost constant as a function of C^o- We notice that the 
theoretical predictions for Bj are typically one or two MeV larger than experiment. This 
is not surprising as it is well known that the inner fission barrier is affected by triaxiallity 
and its height typically decreases by one or two MeV when triaxiallity is included in the 
calculation. Triaxiallity is not included at present because we still do not have access to a 
triaxial code incorporating the BCPM functional but work in this direction is in progress. 
The situation is slightly better in the comparison with the En and S/jvalues. For them, no 
triaxial effects are expected and the agreement is better than for the Bj. 

In table[TT]the theoretical results obtained for different choices of the Eq and i] parameters 
are compared to the known experimental values. The experimental tsF values span a range 
of 27 orders of magnitude for a mass range A = 232 — 286. The theoretical predictions, not 
including triaxial effects and computed with the GCM masses and zero point energies, span 
an even larger range of values and show a large variability depending upon the choices for the 
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Nucleus 


S/(th) 


En (th) 


B n (th) 


Bi(exp) 


-E//(exp) 


^//(exp) 


234 U 


5.87 


1.78 


5.59 


4.80 




5.50 


236 tj 


6.49 


1.90 


6.04 


5.0 


2.75 


5.67 


238 TJ 


6.99 


2.03 


6.54 


6.30 


2.55 


5.50 


238p u 


6.91 


1.85 


5.20 


5.60 


2.4 


5.10 


240p u 


7.43 


2.08 


5.69 


6.05 


2.8 


5.15 


242p u 


7.72 


2.27 


6.30 


5.85 


2.2 


5.05 


244p u 


7.89 


2.47 


6.30 


5.70 




4.85 


240 Cm 


6.8 


1.2 


3.90 




2 




242 Cm 


7.4 


1.7 


4.5 


6.65 


1.9 


5.0 


244 Cm 


8.0 


1.9 


5.0 


6.18 


2.2 


5.10 


246 Cm 


8.4 


2.3 


5.5 


6.0 




4.80 


248 Cm 


8.34 


2.04 


5.47 


5.80 




4.80 


250 Cf 


8.65 


1.25 


4.24 






3.8 


252 Cf 












3.5 



Table I: Fission barrier height parameters Bj (inner) and Bn (outer) as well as excitation energy of 
the fission isomer En. The three parameters are given in MeV. The theoretical values have been 
obtained from the rotational energy corrected HFB potential energy surface. The experimental 



values are taken from [39] for the En and from 



40|| for the B's. 



parameters. Focusing on the "standard" theoretical values 77 = 1.0 and E = 1.0 we observe 
differences of up to 16 orders of magnitude for the lighter nuclei that decrease to differences 
of just a couple of orders of magnitude for the heavier ones. The largest differences are 
observed for nuclei with higher and wider barriers where the impact of parameters like Eq is 
larger. The comparison in isotopic chains indicate that the behavior with neutron number 
is not as bad as the absolute values. The same conclusion can be extracted from the overall 
trend with mass number obtained from the table. Therefore, we conclude that the HFB 
predictions although subject to large uncertainties due to uncontrolled approximations in 
the evaluation of the different parameters can be used to guess with a reasonable precision 
the trends of tsF with mass number. 
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Spontaneous Fission Half-Life Times (s) 






E = 1.0 MeV 


£ = 1-5 MeV 


experimental 




Nuclei 


77=1.10 


77=1.05 


?7=1.00 


77=1.0 


arXiv-1302.7117 


1 


232U 


3.842e+23 


4.213e+29 


2.275e+37 


6.761e+30 


2.5e+21 


2 


234U 


7.211e+23 


1.180e+30 


9.987e+37 


1.621e+31 


4.7e+23 


3 


236U 


8.697e+23 


2.205e+29 


9.759e+34 


1.744e+28 


7.8e+23 


4 


238U 


4.018e+23 


7.371e+28 


7.117e+34 


7.160e+28 


2.6e+23 


5 


240Pu 


3.107e+14 


8.901e+19 


3.567e+26 


1.269e+21 


1.5e+18 


6 


248Cm 


8.684e+6 


1.798e+ll 


3.707e+16 


8.413e+12 


1.3e+14 


7 


250Cf 


2.747e+5 


3.534e+10 


1.074e+17 


2.048e+13 


5.2e+ll 


8 


250Fm 


1.719e-l 


2.937e+3 


5.098e+8 


3.742e+4 


2.6e+7 


9 


252Fm 


1.006e00 


3.971e+4 


4.810e+9 


3.189e+5 


4.0e+9 


10 


254Fm 


4.648e-5 


1.941e-l 


1.181e+3 


5.037e-2 


1.9e+7 


11 


256Fm 


3.327e-6 


1.932e-3 


2.910e+l 


2.447e-2 


1.0e+4 


12 


252No 


1.090e-6 


9.505e-3 


2.782e+2 


8.740e-2 


1.2e+l 


13 


254No 


5.490e-5 


6.689e-l 


3.377e+4 


1.118e+l 


3.0e+4 


14 


256No 


3.853e-5 


1.097e-l 


1.424e+2 


2.071e-2 


l.le+2 


15 


256Rf 


1.355e-5 


1.050e+l 


1.277e+8 


3.337e+4 


6.4e-3 


16 


258Rf 


3.157e-5 


2.879e00 


9.416e+5 


3.892e+2 


1.4e-2 


17 


260Rf 


5.623e-5 


6.264e-l 


6.856e+4 


4.015e+l 


5.1e-2 


18 


258Sg 


1.042e-7 


5.609e-5 


9.134e-2 


1.611e-3 


5.2e-3 


19 


260Sg 


8.999e-7 


6.061e-4 


8.862e-l 


1.517e-2 


7.2e-3 


20 


262Sg 


9.885e-6 


8.554e-3 


2.169e+l 


4.045e-l 


7.0e-3 


21 


264Hs 


3.585e-7 


8.266e-4 


2.735e+l 


5.065e-l 


1.6e-3 


22 


286F1 


2.502e-10 


1.039e-5 


2.625e00 


1.724e-3 


1.3e-l 



Table II: Spontaneous fission half lives for a set of nuclei. Different theoretical predictions corre- 
sponding to different choices of Eq and r] are given. The experimental data is taken from 
text for discussion. 
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D. Neutron rich uranium isotopes 

In the previous section we concluded that the description of fission based on the HFB 
theory is subject to large uncertainties coming from a poor understanding on the way the 
different quantities should be computed. However, we also concluded that the HFB theory 
is giving reasonable well the tsF trends with mass number. Based on this conclusion, we 
have performed calculations in the uranium isotopic chain from the light uranium 226 U up 
to the neutron drip line corresponding to 282 U with the aim of understanding and analyzing 
the trends in spontaneous fission half lives. In Fig[3]the HFB potential energies are depicted 
as a function of Q 2 o for the isotopes of uranium 240 U, 250 U, 260 U, 270 U and 280 U. The energies 
of the heavier isotopes have been shifted by different amounts of energy (0, 55, 100, 135 and 
160 MeV, respectively) to fit all the plots in a single plot. 

We observe the ground state evolution from a quadrupole deformed ground state in 240 U 
with 02 — 0.26 to an spherical one for 270 U (corresponding to N = 178) that is still spherical 
for 280 U. It is also worth mentioned the existence of a second fission isomer in 240 U and 250 U 
that disappears in 260 U, where the first fission isomer is shifted to larger quadrupole moment 
values and zero octupole moment. Two fission isomers reappear in 270 U but both are located 
at a very high excitation energy and different quadrupole deformations than the ones in the 
light uranium isotopes. As a consequence of the higher and wider barriers as the neutron 
number approaches the neutron drip line we expect increasing tgp values as can be observed 
in table IIHI A decrease in the tsF values is observed up to mass number A = 256 but from 
there on a steady increase with mass number is observed. The increase in tgp is very large 
for 262 U that corresponds to a neutron number of 170. From this isotope up to the last one 
considered the values of tsF are very large exceeding in more than 20 order of magnitude the 
typical values of the isotopes close to the stability line. As in previous cases, the tsF values 
have been obtained with the GCM collective mass and not exploring the consequences of 
triaxiallity in the first barrier. 

IV. CONCLUSIONS 

The fission properties of several actinides and super-heavy nuclei have been computed 
with the EDF BCPM recently proposed. The theoretical results for the spontaneous fission 
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Uranium Spontaneous Fission Half-Life Times (s) 





E = 1.0 MeV 


E = 1-5 MeV 


Experimental 


Nucleus 


r7=1.10 


77=1.05 


77=1.00 


r?=1.0 




226U 






1.193e+26 


1.516e+19 




228U 






1.678e+33 


1.622e+26 




230U 






1.724e+35 


3.739e+28 




232U 


3.842e+23 


4.213e+29 


2.275e+37 


6.761e+30 


2.5e+21 


234U 


7.211e+23 


1.180e+30 


9.987e+37 


1.621e+31 


4.7e+23 


236U 


8.697e+23 


2.205e+29 


9.759e+34 


1.744e+28 


7.8e+23 


238U 


4.018e+23 


7.371e+28 


7.117e+34 


7.160e+28 


2.6e+23 


240U 






8.151e+34 


2.059e+29 




242U 






4.425e+35 


1.109e+30 




244U 






3.414e+32 


4.121e+26 




246U 






5.066e+28 


3.068e+21 




248U 






1.266e+26 


3.943e+19 




250U 




2.601e+23 


2.776e+27 


9.389e+19 




252U 






4.862e+29 


1.846e+21 




254U 






1.796e+29 


5.970e+21 




256U 






5.948e+27 


2.473e+20 




258U 






1.786e+29 


3.608e+22 




260U 




1.058e+26 


4.071e+31 


2.577e+24 




262U 






3.416e+43 


1.952e+36 




264U 






4.781e+49 


4.567e+42 




266U 






4.921e+58 


7.239e+51 




268U 






5.448e+68 


7.013e+61 




270U 




4.314e+69 


3.123e+76 


1.317e+68 




272U 






1.126e+83 


1.680e+75 




274U 






7.609e+84 


6.402e76 




276U 






5.363e+98 


7.268e+89 




278U 






1.600e+86 


3.652e+78 




280U 




1.926e+69 


4.209e+74 


1.109e+67 




282U 






3.£|2e+63 


3.351e+56 





Table III: Spontaneous fission half lives for all the uranium isotopes from the proton (A = 226) to 



half lives show a large variability consequence of uncertainties in the evaluation of some 
parameters of the theory and also on the strong dependence of the collective inertia with 
pairing correlations. As a consequence of the large uncertainties in the theoretical results we 
are only able to compare with the experimental data trend with mass number (for instance 
the reduction by 27 orders of magnitude in the spontaneous fission half lives in going from 
A = 232 to A = 286). The theoretical predictions seem to reproduce such trend giving 
us confidence in the relevance of the study of fission barrier properties of neutron rich 
uranium isotopes. There we find that the spontaneous fission half lives remain more or less 
constant up to A=260 where they increase enormously as a consequence of the proximity 
to the neutron number N=170 that is predicted to be a new magic number. On the other 
hand, a comparison of the parameters defining the potential energy surface for fission (inner 
and out barrier heights and fission isomer excitation energies) with the model dependent 
"experimental data" show a rather good agreement that gives us additional confidence on 
the validity of our conclusions. The results obtained clearly show that more attention has to 
be paid to a proper description (including beyond mean field effects) of pairing correlations 
in the configurations relevant to fission. 
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Figure 2: Same as figure Q] but for different pairing strengths. The pairing strengths are given 
in terms of the reference value and a scaling parameter ij taking the values 1.00 (the standard 
calculation), 1.05 and 1.10. 
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Figure 3: HFB energies for some neutron rich uranium isotopes as a function of the quadrupole 
moment. 
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